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Abstract 

We present spectra and mass composition of cosmic rays incoming to the Earth in the energy 
range (0.5— 2)-10^*' eV. As their sources we consider Seyfert galaxies located at distances <40 Mpc, 
following an acceleration model for such moderate-power objects. Mass composition of the particles 
at sources is assumed to be mixed. Generation spectra are described by a function E~^° , where 70 
is an arbitrary parameter. It is shown that the assumptions adopted make it possible to describe 
experimental data provided by HiRes and Pierre Auger Observatory, using different values of 70 . 

1 Introduction 

Since the Greisen-Zatsepin-Kuzmin (GZK) cutoff for protons virith energies >5-10^^eV was pre- 
dicted in 1966 [U [2], issues regarding ultrahigh-energy cosmic rays (UHECRs) have remained one 
of the main topics in high-energy astrophysics. Because of vanishingly small flux of such energetic 
particles, arrays enormous in area should be used over a long period of time. At present a sufficient 
number of experimental events have been collected to provide us with an idea on the behavior of the 
differential energy spectrum of UHECRs. Now the existence of the GZK cutoff is generally considered 
to be proved by measurements of the largest UHECR facilities, HiRes |3] and Pierre Auger Observa- 
tory (Auger) [3]. It is appropriate, however, to recall that the experimental data reported by AGASA 
[5] indicated the opposite effect. 

Although the suppression of the CR flux at ultrahigh energies has been reported by both HiRes and 
Auger, their spectra differ significantly in shape. For instance, while the Hires collaboration clearly 
locates the start point of the GZK cutoff [6j, one of two fits of the Auger spectrum has a smooth 
form ff]. Below we discuss the approximations suggested by both teams. Another, more surprising, 
contradiction between their observations is that the HiRes results demonstrate the proton-dominated 
CR composition at energies from 2-10^^ eV up to 5-10^^ eV [8], whereas Auger points out a transition 
from protons to heavy nuclei, starting with energy of a few times 10^^ eV [9J. This problem was 
discussed, for example, in [TU]. 

To describe experimental data and therefore to determine likely sources of UHECRs, various 
models of initial mass composition and acceleration conditions were considered (see, for instance, 
|llj - [15] ) . Nevertheless the origin of ultrahigh-energy particles remains an open question. In this 
paper we present spectra and mass composition obtained within a model |16| [T7] in which nuclei 
of moderate-power Seyfert galaxies are put forward as UHECR sources. Comparison of calculation 
results with measurements allows us to reveal appropriate acceleration regimes and besides to define 
somehow distances to the most probable sources. It is well known that particles generated above 
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the GZK threshold shift to the region of lower energies in propagating through the space, shaping 
peculiarities in the cosmic ray spectrum, the so-called dip and bump [TH] - 120]. Study of these spectral 
features provides insight into UHECR sources. But a more natural approach for investigation of 
UHECR particles is to exploit data near their maximum acceleration energy, so here we restrict our 
consideration to the region above the GZK cutoff. 

2 Assumptions 

Main features of the above-mentioned model are as follows. Charged particles are accelerated 
to ultrahigh energies at the front of shocks generated in relativistic jets of sources. The jet matter 
resembles one of the accretion disk, therefore the accelerated particles can be both protons and heavy 
ions. Because of lack of information about magnetic field strength within the jets, its value can 
be regarded as an arbitrary parameter specified on the assumption that particles should reach their 
maximum acceleration energies. (Though magnetic field within the jets of some active galactic nuclei 
has been estimated, e.g. [2l| |22].) The expression of the maximum energies for nuclei with charge 

2, derived on basis of this model, has the form: 

^^max = 9.3 • 10^9^2/3 eV. (1) 

The dependence Z^/^ arises due to energy losses (see [Ml HZ]). For protons it is necessary to take into 
account A = Z = 1 so their maximum energy is equal to 

EP,,«4-10^9eV. (2) 

In Section [5] the numerical values in Eq. ([T]) and ^ are varied a bit to obtain a better agreement with 
the experimental data. Such operation is fully justified as at the moment it is hardly possible to set 
these values unambiguously. 

The mass composition of the particles at sources is assumed to be mixed. We consider stable nuclei 
with the highest natural content. Besides protons, these are nuclei of ^He, ^■^C, ^^N, ^^O, ^'^Ne, ^^Na, 
24Mg, 27A1, 28si, 32g^ 40Ar, 40^^^^ 52q^ 56pg_ ^heir abundance ratios are taken from 

The model proposed makes it possible to estimate the maximum acceleration energies -Emax but 
does not specify the shape of the energy spectrum near them. So we adopt a reasonable assumption 
that particle energy spectra at sources obey a broken inverse power law with indices 70 at energies 
E ^ ii^max and 7o+4 otherwise. A similar form was used in [24] to describe roughly the energy spectrum 
of cosmic rays accelerated in supernova shocks. It is necessary to stress that at large values of 70 the 
broken power law should be used only in the vicinity of -Bmax as its extrapolation to substantially 
lower energies may demand too high efficiency from the sources. In the case of a power law spectrum 
the number of nuclei with mass number A and energy in the range of {E,E+dE) is proportional to 
A'^°~^ [25]. If the injection index is high enough, the fraction of heavy nuclei is essentially increased 
at a given energy. The generation spectra of some abundant nuclei at energies of (0.5 — 4) • 10^^ eV are 
presented in Figure [T| with two different values of the index. For 70 = 2.6, nuclei of helium dominate 
among the particles up to ~2.5-102''eV, while nuclei of oxygen and iron prevail alternately, starting 
with £^ « 3- 10^*^ eV. If 70 = 4.6, the pattern changes substantially and iron nuclei begin to dominate 
over the whole energy range of interest. 

The acceleration model adopted was developed for moderate-power Seyfert nuclei. As possible 
sources of UHECRs, this type of active galactic nuclei (AGNs), with redshifts z< 0.009, was identified 
in [26\ I27j . We follow this hypothesis and refer to them as nearby sources. Their distances are set 
according to [28] . 

Based on these assumptions we calculate spectra and mass composition of cosmic rays incoming to 
the Earth at energies from 5 • 10^^ eV to 2 • lO^*^ eV. In this range the main processes by which UHECR 
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Figure 1: Differential energy spectra of some abundant nuclei (p, "^He, ^^O, ^^Fe) at a source, with 70 = 2.6 
(the left panel) and 70 = 4.6 (the right panel). The total spectrum (solid thick curve) is normalized to unity at 

particles lose energy in propagating through the space are 1) disintegration of nuclei due to interactions 
with infrared (IR) and cosmic microwave background (CMB) photons, and 2) pion production by 
protons on CMB radiation. In addition, to estimate the contribution of nuclei to secondary protons, 
we extend the energy range up to Lorentz factor 7Ri2.1-10^^, so photopion production by the nuclei is 
also taken into account. It should be noted that the pair production process is ignored here in order 
to take advantage of the Lorentz factor conservation in photodisintegration and photopion production 
(see Eq. ([3|). Such approximation is reasonable for light nuclei but it may not be reliable for heavy 
species. In fact, the absorption of iron nuclei due to pair production constitutes an essential part 
(~50%) of the total absorption at energies approximately from lO^'^eV to l.T-lO^^eV [29]. But since 
we consider sources of UHECR located mainly at distances below 40 Mpc and the mean free path of 
iron nuclei exceeds 100 Mpc at the energy range mentioned, our assumption is reasonable. 

Computing the mean free paths of particles requires information about IR background radiation 
and cross sections of the above-mentioned processes. Data on the spectral density of IR radiation in 
host galaxies and intergalactic medium are taken from [30] and [31], respectively. Our calculations 
show that interactions of nuclei with galactic background do not lead to any substantial changes, 
thus the results presented below are obtained assuming only the intergalactic absorption. The cross 
sections of different photodisintegration channels for isotopes of nuclei from Li to Fe are computed 
with the TALYS code |32] . For nuclei ^H, ^H, ^He and "^He, we use the total photodisintegration 
cross sections from the Geant4 code [3H]. Each of them includes two reaction channels assumed as 
equiprobable. The photopion production cross sections for ^^Fe, ^He and p are also derived by means 
of Geant4. Corresponding cross sections for nuclei with j4>4 are determined by multiplying that for 
^^Fe by A/ 56 as it follows from the concept of a universal curve |341 135j. In the case of A^3 the cross 
sections are obtained by interpolation between those for ^He and p. A comparison between mean free 
paths calculated for proton and iron nucleus and those from |36j shows that they are in reasonable 
agreement at energies under consideration. 

3 Calculations 

In this section we present computational approaches used to describe propagation of protons and 
heavy nuclei as well as to construct their spectra. 
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3.1 Nuclei 

We consider both stable and radioactive nuclei (with lifetimes of more than several minutes) arising 
due to photodisintegration of particles accelerated at their sources (see Section [2]) . Over a hundred 
isotopes in total are taken into account. 

When an ultrahigh-energy heavy nucleus propagates through the space, its energy per nucleon 
remains constant and its total energy is shared between particles in proportion to their masses after 
disintegration. Thus Lorentz factor 7 is the same for all nuclei produced. This is not correct for 
photopion production because pions take some energy away. But this process becomes dominant 
at very high energies when one needs to estimate the flux of secondary protons. Since the number 
of nuclei with such energies within the sources is exceedingly small and, as indicated below, their 
contribution to the secondary protons can be neglected, we assume that the energy is shared only 
between a nucleon and a fragment in proportion to their masses. Note that only one pion production 
is considered. Two pions production begins to contribute substantially at energies > 550 MeV in the 
mirror system (in which the nucleus is at rest), so this assumption is quite reasonable for interactions 
between nuclei with 7< 1.6-10^^ and CMB photons. Lastly, we presume that protons and neutrons 
are produced equiprobably. 

Our problem can be formulated as a chain of differential equations for the secondary nuclei resulting 
from an initial one. All nonhomogeneous terms of any equation at a given argument are defined if 
previous equations of the chain are already solved. Let N^{x) be the number of nuclei with charge Z 
and mass number ^ at a distance x and A* be the reciprocal of the mean free path for a reaction channel 
i. We consider the following channels: photodisintegration with emitting one or more nucleons, a- 
particle, etc. (every isotope accounts for eight channels of photodisintegration on average), beta-decay 
and pion production. Then change of the nuclei flux with distance is described by the differential 
equation: 

_ ATAr^\ \tot 



dx 



Nt^ix) Af^;,_,^ +E<"^(-) (3) 

n=0,l 

where the reciprocal of the total mean free path takes on form 

At°* = E A' + A'^^+ + (4) 

i 

The expression for A' in the case of a nucleus propagating through isotropic background with 
differential density n{e), is given by [371 138]: 

00 00 

A^ = ^ y de'a, {e') e'J de ^ , (5) 

where e'^^ is the threshold energy in the mirror system and Gi is the cross section of the channel i. 
The solution of Eq. ([s]) can be represented by: 

N^{xn+i) = Ni{xn-i) e-^(tz)(-"+i-"-i) ^ 



/ iV#:j(x') e-^(X.)(— (6) 
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where integrals in the last term are calculated with Simpson's rule. This algorithm was used earlier 
in |39j . As a control, we determine the total number of nucleons, which must be equal to the mass 
number of an initial nucleus, at each integration step Xn- Accuracy of the test is better than 10~^ for 
all distances of interest. 

We should emphasize that the solution procedure becomes more complicated if one includes also 
the energy losses due to pair production as the energy per nucleon does not conserve in this case. 

3.2 Protons 

The dominant process at the energies in question, which affects protons during their propagation 
through intergalactic medium, is photoproduction of pions. For the reason stated above we consider 
production of only one pion. In this case the inelasticity coefficient can be derived analytically, given 
that the initial energy Ep of a proton in the laboratory system and the energy e' of a photon in the 
mirror system are known. 

At given Ep a value of e' varies with both the photon energy e and the proton-photon collision 
angle in the laboratory system. Values of this angle are simulated from isotropic distribution. The 
distribution of e' follows from Eq. ([5]): 



After specifying the photon energy e' , the energy and momentum of a nucleon produced are calcu- 
lated in the center-of-mass system. Further they are used to determine the corresponding quantities 
in the laboratory system and to calculate the inelasticity coefficient. All formulae required can be 
found, for instance, in [40j. 

Thus applying Monte-Carlo simulation we can obtain a function Wp{Ep, E, X)dE which represents 
the probability that a proton with initial energy Ep is localized in the range of {E,E + dE) after 
travelling a distance X, and satisfies the normalization condition: 



The number of trial histories in our calculations is equal to 10^. 

Figure [2] illustrates the probability function at Ep= 2-10^''eV and different distances. The curves 
presented are smooth whereas they have pronounced peaks if mean values of the inelasticity coefficient 
are used. 

3.3 Spectra 

The calculation scheme presented above allows us to construct spectra of both primary (i.e. acceler- 
ated at sources) and secondary (i.e. generated as a result of photodisintegration and photoproduction 
of heavy nuclei) protons as well as spectra of nuclei at different distances from sources. 

First we consider the calculation procedure of the spectra produced by protons. The intensity of 
primary ones with energy E aX a distance X can be derived from 
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where Cp is given by 
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Figure 2: Proton probability function W^{Ep,E,X)dE at different travelled distances X, with Sp = 2-10^°eV 
and rf£^ = 5-10^^ eV. The values of points at energy Ep are indicated on the right scale and give the probability 
that a proton conserves its initial energy after propagating a distance X. 



for energies E' ^E^ax, and 

C^{E') = {E^^/E^^r'^o^E'/E^^)-('ro+^) (11) 

otherwise, with E'min = 5 • lO-*^^ eV. This factor sets the number of particles in the injection spectrum 
at a given energy. 

In the case of secondary protons the so-called source function is required: 

S^{E', x) = (IN^^iE', x)/dx, (12) 

where Np{E',x) yields the amount of the protons with energy E' , produced in travelling a distance x 
by an initial nucleus n with mass number An- Then their intensity is described by 

X oo 

I^{E, X)=Jdxx J Cn{E') S'^{E', x) Wp{E', E,X-x) dE' (13) 

E 

and the total intensity is the sum over all nuclei accelerated: 

I^^\E,X) =J2znAl°-'i;{E,X)/j2znAl^-\ (14) 

n n 

Since the energy of the initial nucleus is An times as great as that of the nucleons produced, the factor 
Cn can be written as 

Cn{E') = {AnE'/Ernin)-^° (15) 

if AnE'^E'^^, and 

Cn{E') = {E^.^/E^^r^" {AnE'/E^^r(^'>+') (16) 

otherwise. 

Nuclei spectra can be constructed in a similar way. It is conventional in cosmic rays physics to 
consider groups of nuclei. In our calculations we also combine them in several groups depending on 
their mass number. These arc the helium group He (^ = 3 — 4, (^hc) =4) and the groups L, M, H, VH of 
light (A = 5-ll, (Al) = 8), mean (A = 12-19, (^m) = 16), heavy (^ = 20-39, (Ah)=32) and very heavy 



6 



(^ = 40—56, (j4vh) =48) nuclei, respectively. Here {A) is the mean mass number. In fact each group at 
a given energy E = 'y{A)mp contains nuclei with energies in the range of jAmp^E ^^Amp, where A 
and^ are the minimal and maximum mass numbers of the group. Nevertheless this procedure is quite 
justified as there exist substantial uncertainties in measuring energy and type of incoming UHECR 
particles. Then the intensity of a group g is defined as 

Ig{E,X) = Y,^nAl'-^C3r,{E)Nn{E,X) /Y,^nA:^'^\ (17) 



where Nn{E, X) is the sum of values N^{X) (calculated according to Eq. ([s])) for the secondary nuclei 
with energy E produced by an initial particle n at a distance X from its source and the index n runs 
over all nuclei contributing to the group. For example, the group VH of very heavy nuclei includes 
particles which could be generated during propagation of ^"^Ar, ^''Ca, ^^Cr and ^^Fe. For initial nucleus 
^^Cr the value Nn consists of all terms between Nfg and iV||. The factor Cn depends on the mean 



mass number {Ag) of the group g and is defined similarly to (15) -(16) with replacement of AnE' by 

AnE'/{Ag). 

The spectra obtained in this way are expressed in relative terms, with the total number of injected 
particles at -E'min= 5-lO^^eV taken as unity. 



4 Spectra and mass composition of UHECRs from a single source 

In this section we present spectra and mass composition of cosmic rays in the range of energies 
from 5-lO^^eV to 2-10^''eV, produced by a single source located at various distances X. Though 
the CR spectrum is modified during propagation of the particles through the space, its shape at the 
observation level strongly depends on the injection index 70. Figures [3]-[5] illustrate the differential 
spectra of protons and nuclei groups, produced by the source with 70 = 2.6, 3.6, 4.6 at X= 5, 40, 100, 
400 Mpc. The maximum acceleration energies of the particles are evaluated from Eq. ([T]) and ([2]). The 
total flux at -Emin is less than unity, which results from transition of more energetic protons to lower 
energies and disintegration of nuclei into fragments during propagation. 

All these pictures have some common features. Firstly, the helium group decays rapidly since the 
mean free path Atot of these particles is small enough. Besides, there is no essential replenishment 
from heavier nuclei. A noticeable flattening in the spectrum of the He group below -E^ax arises from 
the change of dominance between two interaction mechanisms: photodisintegration and photopion 
production, accompanied by a slight increase in Atot- Secondly, although nuclei of the L group (Li, Be, 
B) are absent at the source, they appear in small numbers in propagating of initial nuclei. Thirdly, 
as discussed above, one can neglect the contribution of heavier nuclei to the proton flux. The highest 
proton yield is supported by the He group, as a consequence of its fast decrease. As can be seen, the 
fluxes of all nuclei groups undergo the cutoff at sufficiently large distances from the source. This effect 
should be more pronounced if their absorption due to pair production is taken into account. 

At a given value of the injection index 70 it is easy to estimate the CR composition. Figure [6] 
demonstrates the energy dependence of the mean logarithmic mass of particles composing the afore- 
mentioned total spectra: 

(In^) = Y,h{E)HAg)/Y.^g{E), (18) 

a 9 

where g runs over all nuclei groups, including protons. As can be seen, the larger is 70, the heavier is 
the UHECR composition, especially at energies below 10^0 eV. 



5 Comparison with the HiRes and Auger data 

When calculating spectra, we make a simple assumption that all nearby sources are identical in 
cosmic-ray intensity as well as in composition and acceleration conditions for particles. Therefore the 
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Figure 3: Differential energy spectra of the nuclei groups, produced by a single source at distances X= 5, 40, 
100 and 400 Mpc. The flux of secondary protons generated due to disintegration of the helium group is marked 
by p (He). The total spectrum is presented by grey circles. The number of all particles with energy £^ = 5-10^^ eV 
at the source is normalized to unity. The generation spectra are described by a broken inverse power law with 
index 70 = 2.6. 
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Figure 4: The same as Figure [3j but 70 = 8.6. 
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Figure 5: The same as Figure [3j but 70 = 4.6. 
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Figure 6: Energy dependence of the mean logarithmic mass number of particles composing the total spectra 
from Figures [3]-[5] 



contribution of a single source located at a distance X to the total flux should be taken into account 
with weie: ht X~^. Then the average distance to the nearby sources is close to 9 Mpc. In this case 
the UHECR spectrum is modified only slightly and resembles the result obtained for a single source 
located at a distance of 5 Mpc. 

Figure [7] demonstrates the total spectra produced by the nearby sources, with injection indices 
70 = 2.6, 3.6, and 4.6. The maximum acceleration energies of the particles are evaluated from Eq. ([T]) 
and ([2]). The fits to experimental data, provided by HiRes [6j and Auger f7|, are also presented. All 
data are normalized to unity at energy 5 • 10^^ eV. A comparison with the HiRes broken power law 
fit shows that there is no break in the spectrum at 70 = 2.6, which is the most appropriate for light 
composition. In addition, the spectra at 70 = 2.6 and 3.6 do not reveal sufficient suppression. It is 
evident that in the case of short distances the suppression of the UHECR spectrum can arise rather due 
to acceleration limits at the sources than as a result of interactions of CR particles with intergalactic 
background photons. Hence the spectral break at the GZK threshold Egzk may be attained if we 
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Figure 7: Total energy spectra of particles, produced by nearby sources, with injection indices 70 = 2.6, 3.6, and 
4.6. Values of the maximum acceleration energies are evaluated from Eq. ([I]), ([2|. The HiRes and Auger fits 
are shown for comparison. All data are normalized to unity at 5-10^^ eV. 

take it as the maximum acceleration energy of some nuclei group. The group, of course, should be 
dominant at energies below Eqzk to keep the break in the total spectrum at the observation level. 
According to the HiRes results [6], the GZK cutoff starts from approximately 5.7-lO^^eV and the 
broken power law fit yields indices of roughly 2.8 and 5.4 below and above it, respectively. The Auger 
spectrum also testifies for the CR flux suppression but its fit has a smooth shape [71 HI] . In order to 
reproduce the spectral break and light composition simultaneously in the calculated spectra at lower 
70, we set the maximum acceleration energy for protons equal to E'czif = 5.7- 10^^ eV. The values of 
-E-max for the rest of initial nuclei, calculated from Eq. ([T]), are also increased by a factor of 1.4, however 
the Z^/'^-dependence is retained. Below all results are obtained with the assumption of these new 
maximum acceleration energies. 

5.1 Spectra 

Our calculations show that optimal values of the injection index 70 for representing the HiRes and 
Auger data are 2.2 and 4.3, respectively. In Figure[8]the calculated spectra produced by nearby sources 
are compared with experimental points and their fits. All data are normalized to unity at 5-lO^^eV. 
The energy range is extended down to 10^^'^ eV to show the behavior of the measurements at lower 
energies. According to the left panel, our results are described well enough by a broken power law with 
indices of 2.8 and 5.4 below and above the GZK-cutoff energy at 10^^-^^ eV, as predicted by HiRes. 
This agreement is valid up to the energy ~ 10^'' eV whereas at higher energies the above-mentioned 
feature of the helium group affects the calculated spectrum strongly, making it flatten (see Section |4|. 
In contrast, the right panel demonstrates a good accord between the spectrum with 70= 4.3 and the 
smooth function given by Auger, in all the energy range under study. Computational results presented 
in [7] produce this smooth function within a hypothetical model with a pure iron composition. In our 
model the dominance of heavy nuclei results from a large value of the injection index. 

5.2 Mass composition 

In this subsection we discuss the UHECR composition. At present this topic is one of the most 
attractive problems in high-energy astrophysics. In fact, the largest facilities, HiRes and Pierre Auger 



10 




Figure 8: Differential energy spectra of tire nuclei groups, produced by nearby sources with injection indices 
7o = 2.2 (the left panel) and 4.3 (the right panel). The flux of secondary protons generated by the helium group 
is marked by p (He). The total spectrum is presented by grey circles. The HiRes and Auger experimental points 
along with their fits (grin dashed curves) are shown for comparison on the left and right panels, respectively. 
All data are normalized to unity at 5-10^^ eV. 



Observatory, yield conflicting data. The Auger measurements of the mean shower maximum (Xmax) 
and its root-mean-square fluctuations show that mass composition of extragalactic cosmic rays grows 
heavy with energy [9]. On the contrary, the HiRes analysis of the (Xmax) fluctuations indicates 
the proton dominance [8]. Though both compositions are obtained for particles with energies below 
5-lO^^eV, we extend these results to higher energies, assuming that there is no significant change in 
them. 

Figure [9] shows the calculated energy dependence of the mean logarithmic mass number of particles 
composing the total spectra from Figures [s] and 
retrieved from the mean shower maximum data measured at the highest energy Ec by HiRes ^ 
(Ec ^ 5 ■ 10^'^ eV) and Auger |2j (i^c ~ 3.5 • 10^^ eV), by interpolation between the model predictions 
for protons and iron nuclei. We chose two models which provide extreme results: EPOSvl.99 and 
QGSJETOl (see Figure 3 in p]). The color points are shifted relative to Ec by Alg£' = 0.01 for clarity. 
There is no possibility to describe the experimental data of both facilities in the framework of one 
model. Hopefully, this discrepancy will be overcome after the models used for the comparison undergo 
necessary modifications in response to the LHC data. However, if we compare experimental values of 
(Xjnax)) it is evident that the HiRes composition is lighter than the Auger one at ultrahigh energies. 
According to our calculations, the composition at 70 = 2.2 seems to be a mixture of protons and light 
nuclei such as helium, up to 7-lO^^eV, which is consistent with the HiRes data. At higher energies 
it grows gradually heavier. In contrast, the composition at 70 = 4.3 is compatible with the Auger 
results. Thus one can conclude tentatively that mass composition measurements make it possible to 
distinguish between different values of the injection index. Though, at the moment several models of 
hadron interactions are used and their predictions concerning the shower maximum depth dependence 
on energy differ significantly. 



10 For comparison, the values of (In^) are also 
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Figure 9: Energy dependence of the mean logarithmic mass number of particles composing the total spectra 
from Figures |8] (light grey triangles) and 10 (dark grey circles). The results extracted from the highest energy 
experimental values of (Xmax) provided by HiRes (the left panel) and Auger (the right panel) are shown for 
comparison. The points calculated from the EPOSvl.99 (blue) and QGSJETOl (red) predictions are shifted by 
Algii^ = 0.01 apart from their initial energy. 



6 Discussion 

The necessity of nearby sources of UHECR nuclei has been discussed during the last few years. As 
it was mentioned above, Seyfert galaxies located at distances of <40 Mpc were suggested as probable 
sources of UHECRs j27j . The correlation between the arrival directions of ultrahigh-energy events and 
AGNs within 75 Mpc was found in the Auger data \42\ I43j . Specific quantitative constraints on the 
UHECR source population were obtained in ^44j- For uniformly distributed sources with pure silicon 
or iron composition and intermediate cutoff energies (lO^'^'^— 10^^ eV for iron nuclei), the authors found 
that the lack of nearby sources makes impossible to get a good agreement with the Auger data and 
that the nearest sources have to be within 60 Mpc and 80 Mpc for silicon and iron sources, respectively. 

In our work we considered Seyfert galaxies at distances of < 40 Mpc as sources of UHECRs 
and found that it is possible to describe both HiRes and Auger measurements, varying the injection 
index 79. Moreover, diverting from the specific type of sources, we calculated the energy spectra and 
mass composition of particles produced by sources uniformly distributed within a spherical layer from 
1 Mpc to 40 Mpc, with all other factors being the same. As can be seen from Figures [9] and [T0| there 
is also an agreement with the experimental data provided by the facilities, so the results obtained are 
not too sensitive to the distribution form of the nearby sources. 

It is worthwhile mentioning that our calculations are performed in the framework of the widespread 
notion that the influence of extragalactic magnetic fields (EGMF) on the propagation of UHECR 
particles can be neglected. The same assumption is mandatory for identifying the highest energy 
events recorded by HiRes and Auger with astrophysical objects. In fact, Larmor radius at lO^'^eV 
exceeds 400 Mpc even for iron nuclei in homogeneous magnetic fields below 10^^^ G. However, there 
is no generally accepted opinion concerning EGMF. The influence of random magnetic fields with 
B = 10"^'^ — 10~^G and a coherence scale of 1 Mpc was considered in |44j . The authors found no 
additional constraints on the upper bound imposed on the nearest sources distances. But if the 
magnetic field is significantly stronger (10"*^ — 10^^ G) then the UHECR sources must be located at 
distances of a few Mpc, with the only appropriate candidate being the nearest active radio galaxy 
CenA dS]. 

7 Conclusions 

In this paper we have analyzed the propagation of UHECR particles in the intergalactic medium 
from sources to the Earth and obtained their spectra and mass composition in the energy range 
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Figure 10: Differential energy spectra of the nuclei groups, produced by sources uniformly distributed within 
a layer from 1 Mpc to 40 Mpc, with injection indices 70 — 2.2 (the left panel) and 4.3 (the right panel). The 
description of this figure is the same as that of Figure |8] 



of (0.5 — 2) • 10^" eV, based on the following assumptions: the UHECR sources are considered to be 
Seyfert galaxies located within the radius of w40 Mpc; mass composition of the particles at the sources 
resembles space content; relations between maximum energies of different nuclei are retained according 
to the acceleration model suggested for moderate-power sources; generation spectra obey an inverse 
power law with arbitrary index 70. The higher is the value of 70, the heavier is the composition of 
nuclei accelerated. Over a hundred isotopes arising due to interactions with background photons are 
taken into account. 

Our calculations have shown that the conjecture of nearby sources located at distances of less 
than 40 Mpc is compatible with the HiRes and Auger experimental data. Although these data are 
contradictory, it has been possible to represent them in the framework of a unifying approach with 
varying injection index. 

It is worth noting that our results depend crucially on the model parameters. We hope that 
a larger number of UHECR events, which should be available in the near future, will allow to examine 
generation spectra at sources in more detail. 
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